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The notion of process optimized doners is introduced. Tiiese are quantum doners optimized 
for the average fidelity of their joint output state with respect to a product of multiple originals. 
We design 1 to 2 quantum bit doners using the numerical method for finding completely positive 
maps approximating a nonphysical one optimally, via semidefinite programming. We discuss the 
properties of the so-designed doners as well as their relations to those known to be optimal with 
respect to the clone fidelity. 
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I. INTRODUCTION 

A quantum process describes what can happen to a 
physical system when its state changes. Mathematically 
it is a mapping from the set of the initial states of the sys- 
tem to that of its final states. Quantum mechanics pro- 
vides well-defined limitations to the process to be phys- 
ically realistic. Apart from having to be linear, positive 
and trace preserving, it should be completely positive, 
too. This is an important constraint in the design of any 
quantum information processing apparatus. 

Quantum cloning is a celebrated counterexample of a 
physically realistic process, producing identical copies of 
a physical system in a given unknown quantum state. 
This mapping is not even linear, hence it is unfeasible 
physically. Of course, if something is unfeasible, it can 
still be approximated, as described first in Ref. [i| . Since 
that, the topic of cloning achieved a broad coverage in 
the literature (see Ref. 0] for a review), including the 
calculation of achievable fidelities and several designs of 
particular schemes for cloning. These designs of opti- 
mal doners are based on some intuitive physical ideas, 
which makes them well understandable and gives a hint 
for laboratory interpretation. 

It is also known that the methods of semidefinite pro- 
gramming in operations research makes it possible to 
find completely positive maps ideally approximating non- 
physical ones. In Ref. Q, where this idea was introduced, 
an example is given which is a universal shifter [^-l6|, a 
nonlinear operation on a single qubit. Besides of the 
valuable analytical techniques, Ref 3]. gives a complete 
numerical recipe to design arbitrary operations. 

In the present work we utilize this numerical technique 
in order to analyze the design of 1 to 2 quantum bit don- 
ers. It turns out that partly due to the nature of the 
optimization method, the so designed doner is not op- 
timized with respect to the fidelity of the clones to the 
input state, but the average fidelity of the whole (two- 
qubit output) with respect to a tensor product of two 
copies of the original. This is an approach interesting 



per se, even independently of the actual numerical or an- 
alytical optimization method. We discuss in detail the 
relation of the so-designed doners to those designed for 
optimal clone fidelity. The dependence of the doner op- 
eration on the initial state of an ancilla, as well as the 
effect of restricting the possible input states is discussed. 
This paper is organized as follows. In Sectiorllll we re- 
view the method for finding optimal CP approximations 
of nonphysical maps via semidefinite programming, and 
discuss its application to quantum cloning, introducing 
the notion of process optimized doners. In Section Hill we 
present and discuss our particular results. In Section ITVl 
conclusions are drawn and an outlook is provided. 



II. DESIGN METHOD 

In this Section we first briefly summarize the method 
for designing CP maps via semidefinite programming, 
which optimally approximate an unphysical map. This 
is entirely described in detail in Ref. j3|, we just repeat 
the main ideas here for sake of self-consistency. In the 
second subsection, the application of this method to the 
design of quantum doners is discussed, which leads us to 
the concept of process-optimized doners. 



A. Optimizing CP maps via semidefinite 
programming 

Consider a quantum operation S : B{'Hin) '-^ B{'Hout), 
mapping states in the Hilbert space Hin to those of Hout , 
B denoting the set of density operators on the given 
Hilbert space. In order to be physically realistic, the 
mapping S has to be a linear, completely positive (CP) 
trace preserving one. In some cases, however, one might 
want to at least approximately realize processes which 
are non-CP, or not even linear. Exanxples include the 
quantum shifter discussed in Refs. [J-l6| and quantum 
cloning studied here, etc. 



So we consider an ideal process Sid which is arbitrary, 
and we seek a reahstic (hnear, trace-preserving, CP) map 
iS which approximates iSjd to the highest extent. This lat- 
ter should be quantified somehow. In order to do so, 
consider an input set T C S('Hin), containing the states 
for which we would like our approximate process to be op- 
timal. (This possible restriction might be of some use, it 
enables us, for instance, to consider non-universal quan- 
tum doners in this framework.) We assume that it is 
possible to integrate over T according to a suitable mea- 
sure, this will be denoted by Jj, dT. The quantity to be 
optimized will be 
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(1) 

the fidelity of the output state of the realistic process 
with respect to that of the desired ideal process, averaged 
over the considered input states. We shall term this as 
process fidelity in what follows. In the special case when 
the set {Sid{g)\g G T} contains pure states only, that 
is, we expect the ideal process to map the states of the 
target space to pure states only, the fidelity in Eq. ^ 
simplifies to 
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(2) 



as Sid{g) is a one-dimensional projector. Throughout 
this paper we shall treat the latter case only, as we are 
interested in cloning of pure states, which ideally results 
in products of pure states. 

According to Ref. Q, this optimization can be per- 
formed in the following way. We consider a fixed basis 
on both the input and output Hilbert spaces. S is sought 
for in its Choi-representation, in which it is represented 
by a Hermitian, positive semidefinite operator X acting 
on Hhi (8) Hout, and the relation of the output state to 
that of the input reads 
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where T means ordinary (not Hermitian) transpose in 
the fixed basis, and 1 stands for the identity operator. In 
this representation the process fidelity in Eq. ([T|) of the 
process can be expressed as 



where 



J- = tr (XR) , 



R= I dTg^(»Sid{g). 
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(5) 



The input of the optimization is the matrix R encoding 
all the relevant information on the set T and the process 
S to be approximated. The output will be the Choi ma- 
trix X of the ideal process and the maximum value J-'* of 
the process fidelity. Again note that for Eq. (U) to hold. 



Sid{g) has to be a pure state, that is, a one-dimensional 
projector. 

In the next step we fix two orthonormal bases; a and r, 
in the linear space of the Hermitian matrices over 'Hin and 
■Hout respectively. Assuming dimHin = dimHout — d, 
we have d^ basis elements. We chose ctq and tq to be 
proportional to the identity matrix, whereas the other 
basis elements, indexed with positive integers, should be 
traceless. For d ~ 2 the Pauli matrices, for d = 3 the 
Cell-Mann matrices, while for higher dimensions, gener- 
alized Pauli matrices are a suitable choice. 

As derived in Ref. Q, one can construct the following 
semidefinite program (in dual form, that is, in the form 
of matrix inequalities) to optimize T in Eq. ([2]) : 



maximize p — —ex 

subject to Fq + J2i ^i^i — 0, 



(6) 



where > means positive semidcfiniteness, Fo — ^1, and 
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and the indices 1 < i < d^{d^ — 1) are chosen so that 
the j{i)-s take all the possible values from to d^, while 
the k{i)-s take all the possible values from 1 to d^, and 
the relation of the possible (j, k) pairs and is is bijective. 
The vector x is the one to be found while the constant 
coefficients c in Eq. © are defined as 



q = -tr(i?i^). 



(8) 



Note that these coefficients encode the matrix R of 
Eq. (O, which, on the other hand, encodes all the in- 
formation about the process Sid to be approximated as 
well as on the target set T. Having found the optimum 
p* of the semidefinite program in Eq. ([S]), the optimal 
fidelity is given by 



^*^^+^' 



(9) 



whereas the Choi matrix of the process realizing it reads, 
in terms of the vector x corresponding to the optimum: 
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(10) 



This is the recipe derived in Ref. 3] for finding optimal 
approximate realizations of quantum processes. It can be 
used numerically to design approximate processes. As it 
requires a semidefinite solver which is capable of handling 
Hermitian (complex) matrices, SeDuMi 7] appears to be 
a proper choice. To invoke it for a semidefinite program 
formulated in Eq. ^, it is very convenient to use the 
Matlab package of T. Cubitt [8j, which we have done in 
order to develop our Matlab code producing the results 
described in what follows. 



B. Application for quantum cloning 

Consider the case of quantum cloning. In general we 
are given m copies of a physical system, all in the same 
identical quantum state, say |\E'). We would like to ob- 
tain n > m systems, so that the state of each system is 
closest in fidelity to |^). In the simplest case, discussed 
in what follows, that of the 1 — > 2 qubit doners, we have 
a single qubit in state |^), and we obtain two copies in 
states Qi and g2 so that their fidelity with respect to |^) 
is maximal. If we want a symmetric doner, for which 
both fidelities are maximal, we have a problem charac- 
terized by two objective functions. This is not suitable 
for the recipe described in the previous section. Even if 
we consider the fidelity of the two clones to be equal, it 
is not the kind of problem solved above with semidefinite 
programming. 

What we might consider instead is the following. Take 
the system in state |^), and an ancilla, in an arbitrary 
state, say |0). The target set T shall be the set for 
which we are considering to plan a doner, e.g. for a 
universal qubit doner, the surface of the Bloch-sphere, 
but we may consider restricted input sets. Then carry 
out the procedure so that the ideal process should be 
1^) <8) |0) — )• 1^) (8) 1^), which is obviously unfeasible due 
to the no-cloning theorem. The fidelity J^ considered in 
this optimization is, however, not the same as optimizing 
the fidelity of the clones. 

One may realize that it might be even a different prob- 
lem. Namely, we do not look for a process which produces 
two identical copies which resemble the original to the 
highest extent, but we look for a process which produces 
the product of the two originals to the highest extent. For 
a symmetric 1 — > 2 doner, the fidelity of each clone to 
the original is considered, which will be termed as cloning 
fidelity in what follows. The doners we design here, on 
the other hand, are optimal with respect to the process 
fidelity, hence, we shall call them process-optimized don- 
ers. There are two questions addressed here, at least for 
the simplest, 1 — > 2 qubit case. Are the process optimized 
doners optimal also with respect to cloning fidelity? Can 
a doner be optimal with respect to cloning fidelity while 
being suboptimal with respect to process fidelity? 



III. PROCESS-OPTIMIZED QUANTUM 
CLONER DESIGNS 

In this Section we describe our results regarding 
process-optimized 1 — >■ 2 qubit doners designed the 
above-described method. The numerical results are sum- 
marized in Table m 



A. The universal covariant quantum doner 

First we calculate the properties of the qubit- version 
of the universal covariant quantum doner (UCQC) de- 
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FIG. 1. The quantum logic network for the universal quantum 
doner of Ref. 9J , used here as a benchmark. Input ports are 
to to left, whereas the outputs to the right. 



signed by Braunstein et al. [y| We use this particular 
doner design as a benchmark for those designed by us. 
We use this doner since it is known to be optimal and 
it is a particular circuit so all of its properties can be 
investigated. 

The quantum logic network for the UCQC is depicted 
in Fig. [1] The "circuitry" consists of four controlled NOT 
gates. We do not consider all of its capabilities here, just 
one very particular case: If a qubit in a quantum state 
I 'fin) to be cloned impinges on port 1, while the so-called 
program state is chosen to be 
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|0)|1) + |1)|1)), (11) 



then in the outputs 1 and 2 there will be two identical 
clones of the state to be cloned. (Port 3 will hold a 
state related to the input, it will be omitted). Moreover, 
their density matrix will be the mixture of the state to 
be cloned and a complete mixture. This case (i.e. the 
symmetrical mode of this doner) shall be referred to as 
UCQC in what follows. 

An obvious quantity to consider is the cloning fidelity, 
defined as 
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where g^^^ stands for the state of the first clone. (The 
same can be calculated for the second clone, but as we 
consider symmetric doners here, these will be equal.). 
The UCQC attains the optimal value of 5/6, regardless 
of the state. Another quantity introduced qualitatively 
in Section Hi Bl reads 

^P = ((*i„| ® (*in|) e^'"' (l*in)|*in)) , (13) 

the fidelity of the joint state of the two clones with respect 
to the product of two originals. For the UCQC it can be 
calculated, and it will be 2/3, regardless of the input 
state. 

Further quantities of interest may be the entanglement 
as measured by concurrence of the clones or the von Neu- 
mann entropy 



H{g) = -tiglog^g 



(14) 



of the clones and that of the two clones together. The 
concurrence is calculated according to the well-known 
Wootters formula [lOl |: 

C(g(i")) = max(0, Al - A2 - A3 - A4), (15) 
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UCQC, symmetric mode 


1 ^ 0.83333 


1 Ri 0.66667 


i ^ 0.33333 


0.65002 


0.91830 


universal, opt. for ancilla; |0), used ancilla: 0) 


0.83333 


0.66667 


0.33333 


0.65002 


0.91830 


universal, opt. for ancilla; |l, used ancilla; ^1 


0.83333 


0.66667 


0.33333 


0.65002 


0.91830 


universal, opt. for ancilla; |0), used ancilla; -^1 


0.66667 


0.45833 


0.00000 


0.91830 


1.78434 


universal, opt. for ancilla: ^1, used ancilla: 0) 


0.83333 


0.66667 


0.33333 


0.65002 


0.91830 


equator, opt. for ancilla; |0), used ancilla; |0) 


0.83333 


0.75000 


0.33333 


0.65002 


0.65002 


equator, opt. for ancilla; ^1, used ancilla; ^1 


0.83333 


0.75000 


0.33333 


0.65002 


0.65002 


equator, opt. for ancilla; |0), used ancilla; ^i 


0.66667 


0.50000 


0.00000 


0.91830 


1.70058 


equator, opt. for ancilla: ^1, used ancilla: |0) 


0.83333 


0.75000 


0.33333 


0.65002 


0.65002 



TABLE I. Input-state independent parameters of the designed doners. The most relevant issues are typeset in bold. 



the A-s being the eigenvalues of the matrix 
y y p(i2)^(i^)y (3(12) in desending order, whereas 

g^^^^ = Oy ® ayQ^^^^'^ay (g) ay, ay being the second 
Pauli matrix. Having evaluated these quantities for the 
UCQC, we have listed their values are summarized in 
the first row of Table H] 



B. Universal doners 

First we consider the state to be cloned an arbitrary 
one on the Bloch-sphere: 
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hence, the cloner to be designed is universal. We apply 
the method described in Section|lIl Assume that we have 
an ancilla initially in the state ganc- So the whole two- 
qubit state impingement on the apparatus reads 

ftn(0, (b) = |*in(0, </')>(*in(e, 0)| ® ganc, (17) 

and the output for an ideal cloner would be 

|*out,id(^, </))) = |*in(^, </-)) ® |*in(e, 0)) (18) 

This is to be substituted to the matrix R of Eq. ([5]), 
yielding 
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where we average over the surface of the Bloch-sphere. 



The initial state of the ancilla is also an input for the optimization, providing another input to the problem. We 



considered two possibilities. If the ancilla is in the state |0), the matrix reads 
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In spite of the sparsity of the matrices, we list all their 
elements, since their structure is more visible. For ob- 
vious symmetry reasons, one might consider any other 
pure state instead of |0), the results would be equivalent. 

Carrying out the optimization, we obtain the CP map 
given in Appendix [X] Evaluating the cloning behavior 
we found that in case of both of these doners, all the ex- 
amined parameters are input-state independent, thus we 
have designed a state-independent doner, though there 
was no constraint in the optimization to warrant this. 
The parameters of the designed doners are listed in the 
first four rows of Table U for different scenarios: we have 
applied both designs with both of the considered ancillae. 



As the most important consequence, it appears that for 
a universal doner, the process-optimized one's parame- 
ters are equal to those of the UCQC. Hence, this method 
designs an optimal universal symmetric 1—^2 doner. 
The optimality in terms of process fidelity and that of 
cloning fidelity coincide. In addition, the doner is state- 
independent, though we have not prescribed that. We 
have just optimized for average process fidelity. 

Another consequence is that it is possible to design a 
doner which works for a completely mixed ancilla as well. 
Of course, we seek for any CP maps, their realization may 
require several additional ancillae, hence, one may com- 
bine the doner designed for a pure ancilla to a process 



which replaces complete mixture with a pure ancilla. So 
this consequence is maybe less surprising. Nevertheless, 
as one would expect according to the previous argument, 
the doners designed for the complete mixture as an an- 
cilla perform optimally for the pure input state as well, 
while those designed for the pure ancilla operate as clas- 
sical copiers (cloning fidelity of |), though the clones are 
unentangled. 



It is interesting though to take a look at the Choi 
matrices of the above-mentioned two doners, given in 
Eqs. (IA1[) and (jA2p in the Appendix. Interestingly, the 
rank of the Choi matrix (thus the number of Kraus op- 
erators in the orthogonal Kraus representation) is 4 for 



the completely mixed ancilla, while it is 10 for the pure 
one. This suggests that the doner for the mixed ancilla 
is a "simpler" operation than that for the pure one. 

In summary, in case of an 1 — >■ 2 universal symmetric 
qubit doner, the optimization of process fidelity yields 
a state-independent doner which performs, at least in 
terms of the examined parameters, exactly as the UCQC. 



C. Restricted input doners 

Now we restrict our attention to the optimization for 
the equator of the Bloch-sphere. In this case the matrix 
R in Eq. p^ . the integral becomes a line integral: 
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Note that there is very little difference between this ma- 
trix and that of the universal doner with pure ancilla in 
Eq. (EOD. 

Carrying out the optimization we obtain the doner 
given exactly in Appendix |21 Importantly, for the states 
on the equator this doner is also state-independent. Its 
data are listed in the fifth row of Table HI What is impor- 
tant to note that almost all parameters are equal to those 
of the UCQC, except for the purity of the two clones to- 
gether and the process fidelity. For the equator of the 
Bloch sphere, the process optimized one, being still an 
optimal doner, attains a higher process fidelity. 

It is also worth analyzing the behavior of the doner 
designed for the equator on the rest of the Bloch-sphere. 
For symmetry reasons it is sufficient to investigate a 



I 

meridian of the sphere, that is, the dependence of pa- 
rameters on 9 with, say, = 0. These functions are 
plotted in Fig. [21 where the parameters for the UCQC 
are also plotted for reference. It appears that the param- 
eters reach the values of the UCQC for 6 = ^, that is, 
the equator. As for the behavior of the von Neumann 
entropy, the entropy of each clone is equal to that of the 
system of two clones. In the optimal case this reaches 
the entropy of clones of the UCQC, while at the "poles" 
of the sphere it reaches the value of the joint two-clone 
system of the UCQC. 

We have also carried out the same analysis with the 
completely mixed ancilla as for the case of the universal 
doner. As reflected by the last four lines of Table HI the 
conclusion to be drawn is the same as for the universal 
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FIG. 2. (color online) Quantities of the equatorial doner as a function of the azimuthal angle 9 on the Bloch sphere. Fc is 
the cloning fidelity, Fp is the process fidelity, C is the concurrence of the clones, H is the von Neumann entropy. In all the 
figures, the respective quantity of the UCQC is also plotted, these are the straight lines, the curves represent the quantity for 
the equatorial doner. In case of the entropy the entropy of the clone is always equal to the entropy of the bipartite states of 
the two clones. The lower straight line is the entropy of the clone, while the upper is the entropy of the bipartite state for the 
UCQC. 



case: the doner for pure ancilla becomes a classical copier 
for a completely mixed ancilla. It is possible to design 
a doner for the completely mixture, which works for the 
pure ancilla as well. 

Again it is interesting though to take a look at the 
Choi matrices of the above doners, given in Eqs. (jA3|) 
and (|A4[) . It appears that the spectrum of these matrices 
is the same as that of the universal doners. 

Another question to be addressed in the case of a doner 
for the equator, or more generally, for a main circle of a 
Bloch sphere is that of the dependence on the ancilla. 
To investigate this we have considered a case when the 
target set is the equator rotated by a given angle about 
the X axis. Carrying out this analysis for various angles, 
we have found that the resulting doner has the very pa- 
rameters of the above-detailed equatorial doner: the de- 
signed doner is optimal and anisotropic, independently of 
the angle between the ancilla state and the chosen main 
circle. This is illustrated in terms of cloning fidelity in 
Fig. |3j for a main circle rotated by 7r/4. 




FIG. 3. (color online) The cloning fidelity as a function of 
the spherical coordinates on the Bloch sphere of the input 
states, for a doner designed for cloning the equator of the 
Bloch sphere rotated by 7r/4 about the x axis. The upper flat 
plane represents 5/6, the fidelity of the optimal doner. The 
plotted quantity is dimensionless. 



IV. CONCLUSIONS AND OUTLOOK 

We have introduced the concept of process optimized 
quantum doners: quantum doners that are optimized 
for the fidelity of the state of the whole set of output 
states to that of a product of ideal clones. We stud- 



ied their design for 1 — >■ 2 qubit doners numerically, via 
semidefinite programming. In all the studied cases the 
so-designed doners have been found to be optimal don- 
ers in the usual sense (fidelity of clones to the originals) , 
too. However, in the case of restricted-input doners, the 
so-designed ones provide an output state of two optimal 
clones which are, together, closer in fidelity to a product 



of two ideal clones. 

The results naturally raise the question whether the 
process-optimized doners are better for certain quan- 
tum information tasks, such as, e.g. eavesdropping in 
quantum cryptography, than those designed for optimal 
cloning. Also, it would be interesting to examine the re- 
lation of our doners to other particular results, and gen- 
eralize the consideration to higher dimensions and larger 
number of clones. The framework provided here is obvi- 
ously suitable for such a study. Finally, we believe that 
this framework, if considered analytically in more detail, 
could provide a deeper understanding of the no-cloning 
theorem and quantum cloning. 
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Appendix A: Choi matrices of some doners 

In what follows we write rational numbers as matrix elements. The results of our calculations were their floating 
point counterparts within numerical precision, so we use the rational counterparts for sake of clarity. It would be 
possible to verify analitycally that these are the optima indeed, by calculating the size of the duality gap in the 
optimization. 

a. Universal doner with pure ancilla The nonzero elements of the Choi matrix are: 
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The rank of this matrix is 10. Its eigenvalues are 1.0 (multiplicity 2) and 0.25 (multiplicity 



b. Universal doner with completely mixed ancilla The Choi matrix reads: 
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(A2) 



The rank of this matrix is 4, the only eigenvalue is 1.0. 



c. Equator, pure ancilla The Choi matrix reads: 
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Its rank is 10, with the same eigenvalues and multiplicities as the matrix for the universal doner with pure ancilla. 



d. Equator, completely mixed ancilla The Choi matrix reads: 
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Its rank is 4, with the same eigenvalue and muhiphcities as the universal doner with completely mixed ancilla. 



